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Abstract 

We consider a class of stochastic partial differential equations arising as a model for 
amorphous thin film growth. Using a spectral Galerkin method, we verify the existence 
of stationary mild solutions, although the specific nature of the nonlinearity prevents 
us from showing the uniqueness of the solutions as well as their boundedness (in time). 

1 Introduction 

This paper shows the existence of a stationary solution for a stochastic partial differ- 
ential equation (SPDE), where the solutions may not form a Markov semigroup due to 
the lack of uniqueness. We consider the family of equations 

d t u = -diu + vd 2 x u-dl{d x uf + £,{x,t), veR, (1.1) 

for a real-valued scalar u(t, x) with t > and x 6 [0, L], subject to suitable boundary 
conditions (e.g. periodic or Neumann type). The symbol £ denotes a noise process 
which should be thought of as the generalized derivative of some Wiener process to be 
specified later on. 



Equations of the type (1.1) arise in the growth of thin films (see e.g. [RML+00, 
SP94, BS95]). The function u(t, ■) describes the graph of a surface at time t > 0. 
Usually these equations are equipped with a lot of physical parameters, which we set to 
1 for simplicity. In some models an additional additive nonlinear term (d x u) 2 appears. 
We can treat that case too, but the analysis is more involved without contributing much 
to the general understanding of the situation, so we do not present it. 

It is easy to verify that there exists a value v c < such that if v < v c , the equations 
under consideration present a linear instability. We will therefore distinguish in the 
sequel the stable case (v > v c ) from the unstable case (y < v c ). This instability 
is responsible for the formation of hills, which is frequently seen in experiments (see 
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e.g. [MMS99] and the references therein). On the other hand we have a quadratic 
nonlinearity that compensates this instability. Unfortunately, this nonlinearity makes 
it difficult to derive uniform bounds on the solution. Moreover, it is an open problem 
how to establish bounds in case of a two-dimensional square, which is obviously a more 
realistic model than the one-dimensional case we treat in this paper. This scenario is 
similar to the Kuramoto-Sivashinsky equation, where there are no results for truly two- 
dimensional domains. 

One very helpful tool in the analysis is the conservation of mass: 

Remark 1.1 The quantity M(u) := J Q L u(x) dx decouples from the rest of the equa- 
tion. Therefore, we can assume without loss of generality that M(u(t)) = for a solu- 



tion of (1.1 ). The various Sobolev spaces appearing in the sequel should be thought of 



as the orthogonal complement to the constant function 1 of the usual spaces. 



The local existence of unique solutions to equations of the type (1.1) is standard 



for sufficiently smooth initial conditions. But the existence of global solutions is much 



more complicated, and was shown in [BG02| or [BGR02] using a spectral Galerkin 



method. Nevertheless the question of uniqueness of global solutions is still open, even 
in the deterministic case, as it is out of reach to show enough regularity. Therefore the 
equation does not necessari ly gene rate a Markov semigroup, and the standard theory 



for invariant measures (cf. fc)PZ96 |) does not apply. 

We show in this paper that there exists nevertheless an "invariant measure" for 
( |l.l[ ). To be more precise, we construct a stationary solution {u(t), t 6 R} such that 
the distribution V t := £(«(<))[] of u(t) is constant in time. 

Our concept of solutions is a martingale solution of the corresponding mild for- 
mulation. Hence, we allow a change of the underlying probability space and consider 
solutions not of the SPDE, but of the corresponding variation of constants formula. 
Since we use spectral Galerkin methods, our approach is similar to previous results 



(see [FG95] or [ CG94 ]) for the stochastic Navier-Stokes equation. One of the major 
differences is that we are not able to use the theory of Markov semigroups. Moreover, 
we were not able to get any uniform bound (in t) on the distribution of solutions when 
the driving force £ is a space-time white noise. We are able to establish such a bound 
only for stationary solutions. Therefore, we will construct the stationary process u as a 
limit of the unique invariant stationary solutions to the Galerkin approximations. 

The paper is organized as follows. In Section |[ we introduce the spectral Galerkin 
approximation and present our main result. The next Section presents compactness 
results and the proof of the main result. In Sections ^| and we will give a-priori 
estimates for the solutions. The final Section ^ contains technical results. 
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2 Notation and formulation of the main result 

Define the space L 2 := {/ e L 2 ([0, LJ) : J Q L f(x)dx = 0} with standard L 2 -norm 
We define A as the linear self-adjoint operator in L 2 formally given by 

A := -fl£ + veil , 

where the domain of definition D(A) consists of all functions / G H 4 ([0, L]) satisfying 
M(f) = and boundary conditions given by the equation. We will write H 4 = D(A) 
for short. Moreover, it is well-known that A generates an analytic semigroup {e tA }t>o, 
and we use the fractional powers of A to define the standard fractional Sobolev spaces 
H s for s > with dual spaces H~ s . 

In the sequel we will need spaces of functions on the whole real line with values in 
Sobolev spaces. We recall the definition of the space C(R, H s ), which is given as the set 
of all functions such that for any a < b, the restriction to [a, b] is in C([a, b], H s ). We 
say that /„ converges to / in C(R, H s ) if and only if /„ — ► / in any C([a, b], H s ) with 
a < b equipped with the standard maximum norm. We define the spaces L 2 loc (R, H") 
in an analogous way. 

We write {e k } k& ^ f° r a complete orthonormal set in L 2 of eigenvectors of A and 
denote by II at the orthogonal projector onto the subspa ce of L 2 spanned by e-y, .. . , ex- 
Then the TVth spectral Galerkin approximations un of ( |l.l| ) is given by the solution of 

d t u N = Au N + U N dl{d x u N f + U N W , (2.1) 

where A can be interpreted either as an N x iV-matrix or as a differential operator acting 



on the range of H^. When considering (2.1 ), we will always take initial conditions in 



the range of LTat. In this equation, W is the generalized time-derivative of a two-sided 



cylindrical Wiener process W on L with covariance operator Q. (See [ DPZ92] for the 



definition of a cylindrical Wiener process.) We make the following assumption on Q: 
Assumption 2.1 There exist positive numbers au and a constant C such that 

Qe k = a k e k and \a k \ < C , 

for all k > 0. 

Notice that this assumption covers the case of space-time white noise (a k = 1). The 
assumption that Q and A have the same eigenvectors implies that we restrict ourselves 
to translationally invariant noise, which is also called homogeneous in the physics lit- 



erature. This assumption is crucial to verify technical results like Lemma 6.1 

,N 



Since (2.1 ) is actually a stochastic d ifferent ial e quation in R with locally Lipschitz 



coefficients, it is well-known (see e.g. [ ]Has80| ] or [ |Arn74| ]) that it possesses (locally) a 



strong solution. Standard arguments allow to show the following proposition, the proof 
of which will be given in Section ^ below. 



Proposition 2.2 For every N > 1 and any initial condition in tl^L , , Equation (2.1) 
possesses a unique global strong solution. Furthermore, the law of this solution con- 
verges in variation norm towards a unique invariant measure Qjv which has bounded 
moments of second order. 
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Consider processes {ujv(Q}teR given as stationary solutions of the A^th spectral 
Galerkin approximation corresponding to the invariant measure Qn- Hence, un satis- 
fies the following stochastic ODE. 



d t u N = Au N + Ii N d 2 x (d x u N f + W L 



N 



(2.2) 



where W N is given by W N (t) = J2k=i a ^k w k {t) with the {ui fc } fceN being a family 
of independent two-sided standard Brownian motions defined on the probability space 
underlying W . Since is stationary, we have C(uN(t)) = Qn for any t G R. By 
^at' T] we denote the path measure of {uN(t)}t£[o,T], and by Vn the measure for the 
whole process in path space. 



It is well-known (see e.g. [ DPZ92 1) that, for any pair t > to, the process un satisfies 
(with probability 1) the following variation of constants formula: 



u N (t) 



_ Jt-t )A 



u N (t a ) 







2Jt-s)A 



U N {d x u N (s)) 2 ds 



Jt-s)A 



dW N (s). 



to 



(2.3) 

Again, we consider the differential operators either as operators on the range of IT at or 
as N x A^-matrices. 



As our solutions of (1.1) do not have enough regularity, we will focus on mild 



solutions, which are solutions of such integral equations. Our main result is 



Theorem 2.3 Consider equation (1.1) with periodic or Neumann b.c. in the stable and 
only Neumann b.c. in the unstable case. Then the family of measures {Q.n}n£N given 
by Proposition [2.2| is tight on L 2 . 

Furthermore, for any of its accumulation points Q, there exists a probability space 
(Q,J-,V), a two-sided Q-Wiener process W, and a stationary stochastic process 
{u(t)}teR with u € C(R, H^) n L 2 loc (R, H 1 ), such that C(u(t)) = Qfor every t e R 
and such that 



u(t) = e (t - to)A u(t ) - 



2(t-s)A 



die 



(d x u(s)fds 



(t-«)A 



dW(s) (2.4) 



to 



holds for all t > to, V -almost surely. 



We will not focus on optimal regularity, bu t we co uld slightly improve the regularity 
of u analogous to Corollary 3.2 and 3.3 of [ BG02 |. Moreover, we could prove that 
support of the measure Q is concentrated in a smaller space than L 2 , but we are far 
from getting enough regularity to prove pathwise uniqueness. 



In the stable case, it is easily possible to prove an analog of Theorem 2.3 with 
Dirichlet boundary conditions, but we do not enter into details here. 



3 Proof of the main result 



The main step of the proof of Theorem 2.3 is a bound on the logarithmic moments of 
Qn that does not depend on AT. The main technical difficulty is that Ito's formula can 



not be applied to (1.1) since the covariance of our noise is not necessarily trace class. 
We postpone the proof of Theorem 3 . 1 below to sections |] and |^. 
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Theorem 3.1 Let Qn be the measure on L 2 invariant for the Nth Galerkin approxi- 
mation. Then, there exists a constant C independent of N such that 



f log (I + \\u\\ 2 cl ) Q N (du) < C , 



uniformly in N. 

Using this result, we turn to the 
Proof of Theo rem |2~j| The tightness of the family {Qn} follows immediately from 



Theorem |3.1| and the compact embedding of C 1 into L 2 . We choose any accumulation 
point Q of {Qn} and assume without loss of generality that Qn converges weakly to 
Q in the space of Borel measures on L 2 . Denote by Vn the law of the (u niqu e in law) 



stationary process associated to the invariant measure Qn by Proposition 2.2. 

In order to construct the process u appearing in the statement, we first show that 
the family of measures {Vn} is tight (it turns out that it is so on the space C(R, H~ 3 ) n 
L 2 loc (R, H 1 )), and then verify that the limiting process obtained by the usual Prokhorov- 



Skohorod argument really satisfies the integral equation (2.4). 

To prove the tightness of the family {Vn}, we consider ujy as a solution of (2.3) 
with initial condition un{0) distributed according to Qn- We denote by W% (t) the 
stochastic convolution given by 



W%(t)= [ e A(t - s) dW N (s) , 



and we define VN(t) '■= Ujv(£) — Wj[ (t). The reason is that the stochastic process vn 
exhibits trajectories with much more time-regularity than un- The process vn is then 
pathwise a strong solution of the random PDE given by 

d t v N = Av N - d*Il N (d x v N + d x W%f , v N (0) = u N (0) . (3.1) 

We will need the following technical lemma, the proof of which is postponed to 
sections |] and |[ 

Lemma 3.2 Fix e, T > and assume that there exists R > such that 

P(||«jv(0)|| > R) < e forallNeN. 

Assume furthermore that u(0) is independent of the Wiener increments for positive 
times. Then there exists R > independent of N such that 

P(IMIc(0,tj2) + |M| L 2 ( o,T,ff 2 ) > R) < £ for all N 6 N. 

Using this result, we verify the tightness of {V^' T ^} on the space C(0, T, H~ 3 ) n 
L 2 (0, T, H 1 ) in a similar way as in [BG02, Section 5], so we only briefly sketch the 
main ideas here. 

Given e > 0, we look for a compact set K E such that V N ' (K e ) is bounded 



from below by 1 — e for all N. Combining Theorem 3.1 with Lemma 3.2, there ex 



ists R such that, with probability larger than 1 — s, Vn lies in a ball of radius R of 
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C(0, T, L 2 ) n L 2 (0, T, H 2 ) and ||ujv(0)|| < i?. Furth ermore , W% 



W A in CCO.T.C 1 ) 



P-a.s. Using standard compactness results (e.g. [Gat93 



Proposition 1] or [DPZ92, 



Proposition 8.4]) for the integral operator appearing in (2.3), we can check that the 
above bounds imply the existence of a compact subset K\ of C(0, T, i? -3 ) such that 

G -ft^ with probability larger than 1 — e, Since vn is also bounded in L 2 (0, T, H 2 ) 
with high probability, we obtain by an interpolation theorem (e.g. [ VF88| , Theorem 
IV.4.1]) the existence of a compact set K 2 C L 2 (0, T, H 1 ) such that vn G #| with 
probability larger than 1 — e. 

Hence, {P*"" } N eN is tight on the space C(0, T, #- 3 )nL 2 (0, T, i/ 1 ). By the defini- 
tion of the projection H N , we readily obtain the convergence of W% = H n Wa — * Wa 
in C(0, T, H ), as Wa is already in that space. Combining both arguments, we thus ob- 
tain the tightness of the family {'P$' T] }jveN on the space C(0, T, i?~ 3 ) nL 2 (0, T, ff 1 ). 
Since this holds for arbitrary time intervals, it is straightforward to extend this to the 
whole line, so {V N } is tight on C(R, H~ 3 ) n L 2 loc (R, H 1 ). We call V* one of its limit- 
ing measures and we obtain a subsequence {"Pa^} that converges weakly to V* in the 
abovementioned space. 

Now we can use Skohorod's Theorem to obtain a new probability space (f2, T, P), 
a Q-Wiener process W on that space, stochastic processes Uk with laws Vk = ~PN k 
solving (2.3) with LTat,.^ instead of Wn, as well as a stochastic process u with prob- 



ability distribution C(u) = V* such that u k —> u "P-a.s. in L 2 loc (R, H 1 ) n C(R, H~ 3 ). 
Hence, Ufe(f) — > u(t) in H~ 3 , and additionally we have V t — Q for all t G R by our 
initial choice of a subsequence. 

To show that u is actually stationary, we first remark that Uk — > u in C(R, H~ 3 ). 
Hence, for any choice of (ti, . . . , t m ) G R m we readily obtain in the weak convergence 
of measures on (H~ 3 )" 1 that C((uk(t{), . . . , iik(t m )) —> £((u(ti), . . . , u(t m )). Since 
Uk is stationary, this immediately implies the stationarity of u. 



Using the P-a.s. convergence as in [BG02, Theorem 3. 1], it is technical but straight- 



forward to verify that u is actually a solution of ( 2.4 ) with respect to W. This completes 
the proof of Theorem 2.3 . □ 



4 The stable case 



This section provides the postponed proofs of the previous sections in the case of 
strictly negative A. We will discuss the necessary changes in order to cover the un- 
stable case in Section || below. We start with the 



Proof of Proposition |2.2| . The claim follows from [DPZ96, Has80] if we can show that 
there exists a constant C such that E||itAr(i)|| 2 < C uniformly in t. By Ito's formula, 
we have 



d\\u N \\ 2 = 2(u N , Au N ) dt + 2(u N , d 2 (d x u N ) 2 ) dt + 2{u N , dW N {t)) 



N 

■E 

k=i 



ou. dt 



(4.1) 

Since (un, d 2 (d x UM) 2 ) = and A is a strictly negative definite operator, the claim fol- 
lows after integrating ( 4. 1 ) on both sides, taking expectations, and applying Gronwall's 



formula. Notice that the bound on the second momenta obtained with this procedure 
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diverges with N and it remains an open problem to establish a bound independent of 
N for arbitrary solutions. □ 



To prove Theorem [34] for the stable case, we first verify an L -bound. 

Theorem 4.1 Let be the invariant measure on L 2 for the Nth Galerkin approxi- 
mation. There exists a constant C such that 



L 2 



log(l + |M| 2 ) Q N (du) <C, 



for all N. 



Proof. By (34), the L 2 -norm ||vat(£)|| 2 satisfies 

d t |M| 2 = 2{v Nl Av N ) ~ 2(v N ,d 2 x (d x v N + 8 X W%) 2 



A 

L 

N 



= 2{v N ,Av N ) - 2 / d x (d x v N ) dx-A d x v N d x v N d x W A dx 
Jo Jo 

- 2 f d 2 x v N (d x W%) 2 dx (4.2) 
Jo 

< -ll^wll 2 - v\\d x v N \\ 2 + M\d 2 x v N \\ \\d x v N \\ Wd^U 
+ 2\\d 2 x v N \\\\d x W%\\ 2 

<-^II^A r || 2 -HI^|| 2 + 8|| WAr || 2 ||a^||t ; + 4||a^||t 

Using the Poincare inequality and the fact that we consider only solutions with vanish- 
ing mean, we see that there exists a positive constant a independent of N (but depen- 
ding on L) such that 

a t ||^|| 2 <-a||« Ar || 2 + 8||^|| 2 ||a^|| 4 0O +4||a^||t 

We define now for any interval [s, t] the quantity Wff t ^ by 

W { l t] = j\\\d x W%(r)\\^dr . 
As a consequence, we have the following a-priori estimate on the norm of vn- 



\\v N (t)\\ 2 < e' at+w ^ \\v N (0)\\ 2 + C f e -^-)+<« (l + ||c^( S )||i) 

Jo 



ds 



< e-vt+w^ || Ujv(0 )||2 + Ce w ^(W^ t] + t) . (4.3) 
Since un — i>jv + W A and (0) = we obtain for some e > fixed later on: 

\\u N (t)f < (1 + e)\\v N (t)\\ 2 + C\\W%{t)\\ 2 

< (1 + e) e - at+w ^ \\u N (Q)\\ 2 + Ce w ^(W { % +t) + C\\W%(t)\\ 2 . 



The stable case 



8 



Note that the constants may depend on e. 

The problem at this point is that the exponential moment of the random variable 
WrQ t ] is infinite. We therefore take logarithms on both sides, yielding 



log (1 + \\u N {t)\\ 2 ) < log ((1 + e) e~ at+w ^ (1 + \\u N { 

+ Ce w ^{W { l t] +t)+C\\W%(t)\\ 2 



1 



(4.4) 



Using Lemma 5.1 it is now easy to verify that we can apply Lemma 6.4 (with 
a constant K independent of N) to the right-hand side of (4.4), where we take the 
conditional expectation w.r.t. ujv(O). Hence, 

E((log(l + \\u N (t)\\ 2 )) 2 | u N (0)) < (log (1 + || UAr (0)|| 2 )) 2 (4.5) 
+ 2 (e + log(l + e)- at + EW$ t] ) • log (l + \\u N (0)\\ 2 ) + C , 

for some constant C depending on e and on the parameters of the problem, but not on 
N. 

At this point, we choose first t sufficiently small such that 



This can be done uniformly in TV by Lemma 5.1. Then fix e so small such that 

at 

£ + l0g(l +E) < — . 

Taking expectations on both sides of ( |4.5| ) and using the stationarity of Un(£), we have 

C 



Elog(l + \\u N 



< 



at 



for fixed t sufficiently small, therefore concluding the proof of Theorem 4. 1 



□ 



Remark 4.2 Theorem \. 1 establishes only a u nifor m bound (in t) for stationary so- 
lutions of our problem. Estimates like (\A) or (4.5) are not sufficient to get uniform 
bounds for arbitrary solutions, and this question remains open. 

Let us now turn to the 



Proof of Theorem 3.1 Using (X.2) we obtain after integration 

\d 2 x v N (t)\\ 2 dt <C ( |M*)|| 2 (1 + \\d x W%(fi\\i>)to 



+ c / \\d x w%(t)\\tdt + c\\v N 

Jo 



(4.6) 
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Using Young's inequality and the Sobolev embedding of H 1 into L°°, we have the 
bound 

2 /"Ml«2„,__||2 i wJV||4 , Q 



i^ii^^cn^ivir+ciiwi^ 



This yields 



/ \\u N (t)\\ 2 cl dt<C [ \\v N (t)f(l + \\d x W^(t)\\t )dt 
Jo Jo 



C 



\\d x W%(t)\\jdt + C + C\\u N (0)\\ 2 
( e -<*+<« \\u N (Q)\\ 2 + e<« (Wj£ t] + t)) 



< C 



x (1 + HflxW^WHl,) dt + C / II&WCa (OllJ dt + C + C||ujv 

Jo 

< C7((e<« + 1) |M0)]| 2 + e<« + l) + l) (l + W^) , 

where we used ([O]). Using Theorem [4. 1| and Lemma 5.1, we immediately obtain 

Elog 



f \\u N (t)\\ 2 e idt + lj < C. 



Finally, Jensen's inequality and the stationarity of mat yield 



Elog (11^(0)11^ + 1) = E / log (\\u N (t)\\ 2 cl + 1) dt 



< Elog 



\\u N (t)\\ z cl dt + 1 < C 



concluding the proof of Theorem 3.1 



□ 



We now turn to the proof of Lemma 3.2 This proof will not use the strict negativity 
of A and is thus still valid in the unstable case. Since we need this bound only for a 
fixed time interval [0, T], we can bound the terms in a rather crude way. 

Proof of Lemma p/A Define 



TU 



N 



T - SUp \\d x W£(t)\ 

te[0,T] 



4 
oo 



Using the factorization method and Sobolev embedding it is straightforward to check 
t hat EW^ < C uniformly in N. This result is established completely analogous to 
[ BMPS01 , Lemma 5.1]. Note that the uniformity in N is not trivial, as the family 
{n^jjvgN ls not uniformly bounded as operators on L°° or C°. 

Using this and the assumption on ||«jv(0)||, we see that for every e > there is an 
R > such that 



P(W# < R and ||uat(0)|| <R)>l-e 



(4.7) 
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Combining ( |4.7| ) and ( |4.3[ ), we see that with probability larger than 1 — e one has 

IMi)|| 2 < e w "\\v N (0)\\ 2 + Ce w "(W? +T) 
< e R R 2 + Ce R (R + T) =: R' , 



for any t € [0, T]. Using (4.6) in the same way, we get 

\\d 2 v N (t)\\ 2 dt < C\\v N (0)\\ 2 + C [ \\v N (t)f(WF + l)dt + CTW? 

Jo 



o 



< C(R + TR'(R + 1) + TR) . 



with probability larger than 1 — e, thus concluding the proof of Lemma 3.2 



□ 



5 The unstable case 

This section deals with the case where the operator A is no longer strictly negative 



definite. In order to treat this case, we make use of a trick that was used in [ NST85 



CEES93Q to get bounds on the deterministic Kuramoto-Sivashinsky equation. It turns 
out that the present model is sufficiently close to that equation to make that trick go 
through. Nevertheless, we can only treat Neumann boundary conditions (which is the 
same as considering the restriction on [0, L] of functions that are even and periodic 
with period 2L). In a similar way we can treat also Dirichlet boundary conditions, but 
periodic boundary conditions are still open. 

Most of the proofs are similar to the previous section, so we will only state the main 
differences. Instead of defining vn as previously, we define vn by 



v N (t) = u N (t) - W%(t) - <$>n , 



(5.1) 



where $jv = IIjv$ for some function $ to be chosen later and is the stochastic 
convolution defined in the previous section. The stochastic process vn then satisfies 
the following random PDE: 



d t v N = Av N + A$jv - n N di(d x v N + d x W% + d x Il N $) 



(5.2) 



We can rewrite this as 

d t v N = Av N - d 2 x (d x v N + d x W%f - dl{dJ\ N Vf + A$ - 2fl£(Q B $ d x W A ) , 
where the operator A is defined as 

Av = Av - 2dl{d x $d x v) . 

Using exactly the same technique as in the previous section, we see that in order to 
get uniform bounds on the Galerkin approximations of v, it suffices to find a smooth 
function $ such that 



(v, Av) < -c\\d 2 x vf 



(5.3) 



for some constant c > 0. Using this function, it is easy to verify that the assertions of 



Proposition 2.2 and Theorems 3.1 and \. 1 hold in the unstable case, too. The only major 
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changes appear in the values of the constants, which do now depend on the choice of 
<£>. We will therefore not go through the proofs of these assertions for the unstable case, 
but we will sketch how to find a function $ such that A is strictly negative definite. 
Integrating by parts, we see that the bilinear form (|5.3[) can be written as 



(v,Av) = -]]d>|| 2 - u\\d x v\\ 2 - (d x v,(d 2 $)d x v) , (5.4) 

where v is negative. The problem is therefore reduced to finding a smooth periodic 
function $ such that the Schrodinger operator 

W» = -§a£ + fi£$(a;), 
with Dirichlet boundary conditions satisfies (u, H<s,u) > | v\ \\u\\ 2 for all functions u in 



its domain. The idea appearing in [NST85] is to choose $ such that, away from the 
boundary, d x Q is for all practical purposes constant and sufficiently large (say equal 



to about 2\v\). The problem is that, in order for (5.2) to hold, $ has to belong to 
D(A) and must therefore satisfy the same boundary conditions as u. As a consequence 
d 2 & must satisfy d 2 <fr(s)dx — 0, which is of course impossible for a constant 



(non-zero) function. Looking at (5.4), we notice that d 2 &(x) — 2\is\(l — S(x)) would 
formally fit our needs, since the delta-peak is integrated against d x v, which vanishes 
at the boundaries, due to the Neumann conditions. The function $ obtained this way 
does of course not belong to D(A), so we look for an approximation of it which is more 
regular. 

Since $ satisfies Neumann boundary conditions, it is natural to write it as 



—J . (5.5) 

n=l 

(The sum starts at 1 because we are interested only in functions with vanishing mean.) 
If we choose ip n = 2ri~ 2 , we see that d 2 & is given by d 2 Q(x) = — 5(x)), which 

is what we would like to approximate. In order to get a regular function, we define 

2 J 2 for n < 2n*, 
^ : =^« = (o forn>2n„ (5 ' 6) 

where n* is some (sufficiently large) constant to be fixed later on. With this definition, 
we have: 

Proposition 5.1 For every L, C > 0, there exis ts a value n* > such that the quadra- 
tic form with $ defined as in f l5.5| ) and ( |5.6| ), satisfies 

L 

u(x)(H^u)(x)dx > C\\u\\ 2 , 
for every u in the domain ofTL^,. 



Remark 5.2 Notice that the function $ defined by (5.5) and (5.6) is actually analytic, 
so the expressions appearing in (5.2) and containing $ can all be bounded uniformly 
in N (not in n* of course, but n* is chosen independently of N). 
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Remark 5.3 As in [ CEES93 ] we could choose a slow decay of ip n for n > 2n* to 



optimize the L-dependence of our bound, but we neglected this for simplicity. 



Proof. Applying the arguments of [ |CEES93| , Prop. 2. 1], we see that it suffices to show 
that the quantity 



£ \4> k+ rn-J k - m \ \ w|(h £i 



k>m>0 



EkE n 



an 



can be made arbitrarily small by choosing rt* sufficiently large. The only non-vanishing 
terms of this sum are those where < k — m < 2n* and k+m > 2n*. We can estimate 
these terms by 



^ 1 2n*+m _^ 



. oo 2n*+m ^ 



a* ' — ' ' — ' m 2 k 2 

m=n» k=m 



a 2 Z-^ Z-^i rrflk 2 a 2 ^— ' ^ 

m— 1 fe-2n*-m m= 

"V— ■ 8 y- 1 



m— 1 " m=n* 



4tt 2 



In both sums, we used the fact that k is larger than n* and that there are less than 2n* 
terms in the inner sum. Thus, Y can clearly be made arbitrarily small by choosing 
n* sufficiently large. This proves Proposition 5.1 and concludes our exposition of the 
unstable case. □ 



6 Technical estimates 

In this section, we prove the two technical estimates required for the proof of Theo- 
rem \. \ above. 



Lemma 6.1 There exists a constant C independent of N such that 

nd*w%(t)\\io < ct 1 ' 8 

for all t < 1. 



Remark 6.2 The power 1 / 8 in the above lemma is not optimal but it is sufficient for 
our needs. All we need is E||9 X W^ = o(l) uniformly in N. 



Remark 6.3 The constant in the above lemma depends only on the coefficients of the 
problem and the bound on the oifc. It is possible to allow for slowly growing otk, using 
the Sobolev embedding of L°° into the fractional Sobolev space W 8,p for sp > 1. 

Proof. For / g L°°([0, L]) with vanishing mean, denote by {/fcjfcgN its Fourier coef- 
ficients. Since the eigenfunctions of A are uniformly bounded in L°°, we have the 
following estimate on ||/||oo: 

OO OO OO j,^ 

ll/l|oo<c£| /fe |<c(]>>|Vv fe | 4 / 3 ) (£l fc r 3/2 ) <c\\ Kf \u, 

fc=l k=l fe=l 
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where K is the operator that acts on Fourier coefficients as (Kf) k — |fc| 3 ^ 8 /fe- Here 
we used the usual isometry between W and £ q forp -1 + q^ 1 = 1. 

Denote by \ k the eigenvalues corresponding to the eigenfunctions e k . By the defi- 
nition of A, there exist constants a such that 

cik 4 < -A* < c 2 k 4 , \{Kd x e k ){x)\ 2 < c z k 14 ' 4 . 

With these notations, [Kd x Wff) (t, x) (with fixed values of x and t) are centered Gaus- 
sian random variables given by 

{Kd x W%){t,x) = Y^a k {Kd x e k ){x) / e~ x ^- s > dw k (s) , 

with independent Wiener processes w k (t). The variance of (if d x W^)(t, x) is thus 
bounded by 

N ,t 

E\(Kd x W%)(t,x)\ 2 = Y / a 2 k \(Kd x e k )(x)\ 2 / e 2X ^ ds 
k=i Jo 

<^Ew 11/4 / e- clfc4 ^<£^(i-e- C2feV )- 
fe=i " fe=i 

We now take some fc* to fixed later and split the sum into two parts: 

fe»-i 



E\(Kd x W A )(x,t)\ 2 <Cj2 WH + CY, T^< Ck * 5/4t +-Tj4 



Choosing fc* w i~ 1/4 , we have the estimate E| (Kd x W%)(x, t)\ 2 < Ct 1 / 16 . Since the 
random variables (K d x W A )(x, t) are Gaussian for all fixed values (x, t), we have 

E||Wf < GE||A-a x W^||t = [ L E(Kd x W A )(x,t) 4 dx < Ct 1 ^ . 

Jo 



□ 



Lemma 6.4 For every e > one/ if > //zere erafs a constant C depending only on e 
and K, such that for any two random variables W\ and W2 with E ( | W\ \ 2 + | W2 \ 2 ) < 
K we obtain 

E(log (xe Wl +e W2 )Y < (\og(x)) 2 + 2(e + EW^logx + C , 
for every x > 1. 

Proof. Expanding the square, we get 



E(log (rre^ 1 + e^ 2 )) = (log a;) 2 + 2 log a: (EWi + Elog (l + e^-^ 
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+ E(Vi + log (1 + e w *- Wl /x 
Since we assumed x > 1, we have 

E(Wi + log (1 + e W2 - Wl /x)Y < 2(EW? + E(W 2 - T^i) 2 + l) 

< 6K + 2 . 



It thus suffices to show that there exists xo > depending only on K and on £ such 
that 

W a -Wi 

E* :=Elog(l 

x 

for x larger than x . 

To verify (5. 1 ) consider arbitrary e > 0. We define y — \Ws — Wi \ and denote the 
probability distribution on R + of y by P. Now choose yo > 1 large enough such that 

poo 

yP(dy) < e , / P(dy) < e . 

yo J yo 

We choose yo — 1 + (Ey 2 )/e. Now define 



e l+2K/e e y 

Xq = > — 



We thus have 



E x = log(l + — )P(dy)+ / log(l + — )P(dy) 

JO V x ' Jyo K X ' 

/ e yo N z- 00 
<log(l + — )+/ log (1 + ee v ~ y °) P(dy) 
x Jyo 

<e+ \og((\ + e)e y - ya )P{dy) 

J Vn 



< e + £log(l + e)+ / (y-y )P(rf2/)<2£ + £log(l + £). 

Jyo 

The claim follows immediately. □ 
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